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Abstract 

A perfectly reflecting (Dirichlet) boundary condition at the edge of an impen- 
etrable magnetic-flux-carrying tube of nonzero transverse size is imposed on the 
charged massive scalar matter field which is quantized outside the tube. We show 
■^j - \ that the vacuum polarization effects give rise to a macroscopic force acting at the 

I increase of the tube radius under the steady magnetic flux. 

O 

o 1 Introduction 

(N' 

Polarization of the vacuum of quantized matter fields under the influence of boundary 
conditions was studied intensively over more than six decades since Casimir [1] predicted 
^ ■ a force between grounded metal plates: the prediction was that the induced vacuum 

energy in bounded spaces gave rise to a macroscopic force between bounding surfaces, 
see reviews in Refs. [2] and [3]. The Casimir force between grounded metal plates has 
now been measured quite accurately and agrees with the theoretical predictions, see, e.g. 
Refs. [4] and [5], as well as other publications cited in Refs. [2] and [3]. 

In the present paper we consider the vacuum energy which is induced by boundary 
conditions in space that is not bounded but, instead, is not simply connected, being an 
exterior to a straight infinitely long tube. This setup is inspired by the famous Aharonov- 
Bohm effect [6], and we are interested in polarization of the vacuum which is due to 
imposing a boundary condition at the edge of the tube carrying magnetic flux lines inside 
itself; this may be denoted as the Casimir- Aharonov-Bohm effect (see also [7]). 
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The temporal component of the energy-momentum tensor for the quantized charged 
scalar field *f>(x) in flat space-time is given by expression 

Too(x) = \ \ + - \ + V 2 (1) 

where V is the covariant spatial derivative involving both affine and bundle connections 
and the field operator in the case of a static background takes form 




A 



a\ and a\ (b\ and b\) are the scalar particle (antiparticle) creation and destruction op- 
erators satisfying commutation relations; wave functions ip\{x) form a complete set of 
solutions to the stationary Klein-Gordon equation 

(-V 2 + m 2 )^(x)=^(x), (3) 

m is the mass of the scalar particle; A is the set of parameters (quantum numbers) spec- 
ifying the state; E\ = > is the energy of the state; symbol ^£ denotes summation 

A 

over discrete and integration (with a certain measure) over continuous values of A. 

As is known for a long time [8, 9, 10], the energy- momentum tensor depends on 
parameter £ which couples \& to the scalar curvature of space-time even in the case of 
the vanishing curvature, see (1); conformal invariance is achieved in the limit of vanishing 
mass (m = 0) at £ = (d — l)(4d) -1 , where d is the spatial dimension. Consequently, the 
density of the induced vacuum energy which is given formally by expression 

e = (vac|T 00 (a;)|vac) = ^>aV'1(x) V'a(x) - (£ - 1/4) V 2 ^X^(x) ^ A (x) (4) 

A A 

depends on £ as well. This poses a question: whether physically measurable effects (i.e. 
the Casimir force) can be dependent on £? 

In the present paper we are considering a static background in the form of the cylindri- 
cally symmetric magnetic flux tube of finite transverse size, hence the covariant derivative 
is V = d — ieV with the vector potential possessing only one nonvanishing component 
given by 

V v = $/2tt (5) 

outside the tube; here $ is the value of the magnetic flux and (p is the angle in polar (r, tp) 
coordinates on a plane which is transverse to the tube. The Dirichlet boundary condition 
at the edge of the tube (r = vq) is imposed on the scalar field: 

V>aU o = 0, (6) 

i.e. quantum matter is assumed to be perfectly reflected from the thence impenetrable 
flux tube. 

As we shall see, the vacuum energy which is induced outside the flux tube gives rise to 
a macroscopic force acting at the increase of the tube radius under the steady magnetic 
flux. Although the induced vacuum energy density is ^-dependent, the Casimir force will 
be shown to be independent of £. 
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2 Vacuum energy density 



The solution to (3) outside the magnetic flux tube can be obtained in terms of the cylin- 
drical functions. The formal expression (4) for the vacuum energy density has to be 
renormalized by subtracting the contribution corresponding to the zero flux. Restricting 
ourselves to a plane which is orthogonal to the tube, we obtain (for details see [11]): 



dk k (k 2 + m 2 ) ^ 2 [S(kr, kr ) — S(kr, kr )\^ =0 ] — 




-(£-1/4) A J dkk(k 2 + m 2 )" 1/2 [S{kr, kr ) - S{kr, kr )\ 



<S>=0 



(7) 



where 

S(kr, kr ) 



[Y\ n _ e ^/2n\{kr Q )J\n-e<5>/27T\{kr) - J| n - e $/ 27 r| (A^o) Y\n-e^/2n\ (kr 



Y 



n-e*/2 



.|(^0) + ^| 2 n _ e$/2 ^|(^ ) 



Z is the set of integer numbers, J M (u) and Y^(u) are the Bessel functions of order \x of the 
first and second kinds, and A = d 2 + r~ l d r is the radial part of the Laplacian operator 
on the plane. 

Owing to the infinite range of summation, the last expression is periodic in flux $ 
with a period equal to 27re -1 , i.e. the London flux quantum (we use units c = K — 1). 
Our further analysis concerns the case of $ = (2n + l)7re -1 when each of the integrals in 
(7) is the most distinct from zero. Introducing function [12] 



G(kr,kr ) = S(kr, kr )\^ = 
we rewrite (7) in the dimensionless form 



S(kr, fcr )|$= 



o- 



re. 



a + (mr , mr) — (£ — l/4)r A 



3 «_(mr , mr) 



where 



a±(mro, mr) 



1 

2^ 



21 ±!/2 



G(z, Az), 



(9) 



(10) 



'111 



and A = r /r (A G [0, 1]). 

We follow the technique of numerical calculations developed in [11, 12] with some mod- 
ifications, notably we perform direct integration over consecutive periods of the G(z, Xz) 
function using the Euler-Maclaurin integration formula [13]. This results in a sufficient 
decrease of the computation time. 

Thereafter we calculate a + and a_ functions for the case of mr = 10 -3 at a set of 
different distances from the axis of the tube. This allows us to obtain coefficients of the 
interpolation function which is approximated in the form 



a±{x ,x) = [±e~ 2x x lTl/2 ] 



X — Xo 

X 



Pt{x 



ar 



Qf{x 2 ) 



x v 



X > Xq, 



(12) 
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Figure 1: Behavior of the a + (x ,x) (positive) and the a_(x ,x) (negative) functions for 
the case of Xq = 10" 3 . The variable x (x > Xq) is along the abscissa axis. 

where x = mr, Xq = mr and P^(y), Q^iv) — are polynomials in y of the n-th order with 
the Xo-dependent coefficients. First factor in square bracket in (12) describes the large 
distance behavior in the case of the zero-radius tube (singular thread), second factor in 
square bracket is an asymptotics at small distances from the edge of the tube, and the last 
factor in square bracket is the intermediate part. Since the flux tube is impenetrable, the 
a± functions vanish at x < Xq. Behavior of the dimensionless a± functions is presented 
on Fig.l. 

For the a + function we estimate the relative error of the obtained result as 0.1%. It 
should be noted that nearly 95 % of the integral value is obtained by direct calculation 
and only nearly 5% is the contribution from the interpolation. The integration in the case 
of the a_ function is performed more quickly and with a higher accuracy, as compared 
to the case of the a + function, because the former tends to zero more rapidly at large 
distances. In this case the contribution from the interpolation can be estimated as 10 _3 % 
from the total value. 

We define function [11] 

- / \ 3 A ot-(x ,x) da-(x ,x) 2 d 2 a-(x ,x) 

a-(x ,x) = r o A = ct-(x ,x)-x - Vx — (13) 

r ox ox 

and present its behavior on Fig.2. 




Figure 2: Behavior of the a_(xo,x) function for the case of Xq = 10 -3 . The region in a 
rectangle on the left figure is seen in the scaled-up form on the right figure. The variable 
x (x > xq) is along the abscissa axis. 




Figure 3: The dimensionless vacuum energy density r 3 e ren (xo, x) at different values of the 
coupling to the space-time curvature scalar for the case of xq = 10~ 3 . The region in a 
rectangle on the left figure is seen in the scaled-up form on the right figure. The variable 
x (x > x ) is along the abscissa axis. 

We construct the dimensionless vacuum energy density at different values of the cou- 
pling to the space-time curvature scalar (£): 

r 3 e ren = a+(x Q ,x) - (£ - l/A)&-(x ,x) (14) 

and present its behavior in the case of mro = 10~ 3 on Fig. 3. 

The behavior of the induced vacuum energy density as the radius of the tube tends 
to zero is of primary interest. To do a numerical calculation at xq < 10~ 3 needs a long 
computational time and is a rather complicated task. Nevertheless, we can make some 
general conclusions regarding the case of small xq. 

It seems plausible that this case with the decrease of the tube radius becomes more 
similar to the case of the tube of zero radius (singular thread) see, e.g., Fig. 4. However 
there are some peculiarities in the behavior in the vicinity of the tube. To discuss them, 
let us first recall the exact expressions corresponding to the case of the singular magnetic 
thread (see [14]): 

a+(0,x) = — 2 | - - 2xK Q {2x) - K x {2x) + ~ J ^ J — (15) 

- 7rx[K (2x)L 1 ( y 2x) + K 1 (2x)L (2x)}\ , (16) 

a_(0, x) = 4 {\ ~ 2xK (2x) - K^2x) - nx [if (2x)L 1 (2x) + K^L^x)]} , (17) 
d_(0, x) = -4[2xA' (2x) + K x (2x)], (18) 

where K v {u) and L u (u) are the Macdonald and the modified Struve functions of order v. 
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Figure 4: The constituents of the dimensionless vacuum energy density: a) a+, b) a_, 
c) a_ for the case of Xq = 10 -2 (dotted line), 10 -3 (dashed line). The region in a 
rectangle on the c-figure is seen in the scaled-up form on the d-figure. The behavior of 
the corresponding functions for the case of a singular magnetic thread is presented by 
solid line. The variable x (x > 0) is along the abscissa axis. 



Consequently, in the vicinity of a thread one gets 

I _ 3-£ 2 

a+(0,x)= ^ , £<1 (19) 
a- 0, x = ^— ^ — + C(x 3 ), x < 1, (20) 



&- (0, x) = - + 1 + ^ + 2 ln x ^ + o (x 3 ) , 



where 7 is the Euler constant. Using the latter relations, we get asymptotics of the 
renormalized vacuum energy density at small distances from the singular magnetic thread 

A* * - £ " f« " V) (~ + 1 + 2 L+ 2hj V l + 0(4 , « 1. (22) 



12tt 2 4tt 2 V V V 2tt 2 2tt 

In contrast to (19) and (20), the a±(xo,x) functions in the case of nonzero radius are 
vanishing quadratically in the vicinity of the tube, see [11], 

a ± (x ,x)\ x _^ X0 ~O[(x-x ) 2 }. (23) 

To be more precise, we assume the asymptotics in the form, cf. (12), 

a±(x ,x) = ±— — jp—f±(x ,x), (24) 
6 
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Figure 5: The expected behavior of the constituents of the dimensionless vacuum energy 
density at small distances from the tube: a) a+, b) a_, c) a_ for the case of Xq = 10~ 9 
(solid line). The behavior of the corresponding functions for the case of a singular magnetic 
thread is presented by dashed line. The variable x (x > xq) is along the abscissa axis. 



then one gets 



d 2 



a-(x ,x) = -(x — x ) —f~(x ,x) + 



! _ 6 ££ + 5 f) 4/_ 

X X z J OX 



[X Q ,X - 



X X' 



l-8^ + 9^| )f ( x ,x), (25) 



with a„(xo,2o) = ~ 2 f-( x o, x ). 

The f±(xo,x) functions are adjusted as 

f+(0,x)= w , 



„, . l-7ra; + (3-27-21nx)x 2 



2tt 2 



consequently, one gets 



(26) 
(27) 



xo -> 27T 2 

a: — >■ xq 



1 l + 27 + 21nx o 4-7rx 
H — z H — x + 



2vr 2 



7T 2 X 



+ 



-9 + Atxx - 7x 2 + 27X 2 + 2x 2 lnx 
2ir 2 x 2 



xl (28) 



The asymptotical behavior of the a± and a_ functions with the use of (24) - (28) is 
presented on Fig. 5 for the case of a sufficiently small value of x . As one can see, this 
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behavior is quite similar to the behavior for the case of Xq = 10 3 , compare with Fig. 2 
and Fig. 4. It should be noted that the f±{xo,x) functions depend strongly on xq. 



3 Total vacuum energy and the Casimir force 

The total vacuum energy which is induced outside the magnetic flux tube of finite radius 
is 



2tt oo 
r-o 

In view of the relation 



oo oo 

a+(x ,x) dx _ft_l\ f a-(x ,x) 
x 2 \ 4 / / x 2 



dx 



■i-t) 



(29) 



a-{xo, x) 



x- 



dx 



ct-(xo,x) da-(xo,x) 



x 



dx 



(30) 



•I'D 



which follows from (13), and relations (24) and (27), we conclude that the vacuum energy 
is independent of the coupling to the space-time curvature scalar (£): 



E = 2nml a+M dx. 



ar 



(31) 



This is in contrast to the case of the singular magnetic thread, when the total induced 
vacuum energy is divergent and £- dependent (see [14]): 



2-7T OO 

E*" = J dtp J e^rdr r.4mU-^\j J. (32) 



It is curious that the vacuum energy in this case is finite at £ = 1/12, being equal to 



o 

- Tix [K (2x)L 1 (2x) + K 1 (2x)L (2x)]^xdx = -0.01989 x 2nm. (33) 

Although vacuum energy E (31) is finite, its absolute value grows infinitely as xq tends 
to zero (see (24) and (26)): 

(34) 



E \xo-*0 ~ 



in 



which is in accordance with the divergence of the vacuum energy in the case of the 
singular magnetic thread. To be more precise, relation (30) fails to yield zero in the case 
xq = 0, and, therefore, the divergence of the vacuum energy in the latter case becomes 
^-dependent. 
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Figure 6: The total vacuum energy as a function of the tube radius in the range 10 -3 < 
Xo < 3/2. The variable x$ is along the abscissa axis, the value of hx-^— is along the 
ordinate axis. Solid line interpolates the dots that have been calculated. 



We present the values of vacuum energy E (31) for several values of the tube radius 
in the Table: 



Xo 


3/2 


1 


1/2 


io- 1 


10~ 2 


IO" 3 


E/(2nm) 


5.013 • 10~ 12 


6.944 • 10~ 10 


2.068 ■ 10~ 7 


1.65 • 10~ 4 


0.0106 


0.1486 



These results are also given on Fig. 6 in logarithmic scale, where the dots corresponding 
to the data in the Table are joined with the help of the interpolation function, 

E 

r){xo) = In , (35) 



2irm ' 



which can be taken in the form 



77(2:0) =a + ^2 



1 



lnx - 



(36) 



where bi and dj are the positive adjustable constants. 

To change the radius of the magnetic flux tube one has to apply a work that is equal 
to the change of the total vacuum energy which is induced outside the tube. In the case 
of the infinitely small change of the radius one has 



AE = 2nproAro, 



(37) 



where p can be interpreted as the vacuum pressure which acts from the outside to the 
inside of the tube 



1 dE 



3 »?0o) 



m 



dr](x ) 



27ir dr xq dxo 

This results in the Casimir force acting from the inside to the outside of the tube 



F(x ) = -27rr p{xo) = -2ixm 2 e vixo) - 



dxo 



(3* 



(39) 



The behavior of the Casimir force is presented on Fig. 7. 

As one can see, the Casimir force tends to increase the radius of the tube and to min- 
imize the induced vacuum energy of the quantized scalar field. Certainly, our conclusion 
is obtained under the assumption that the magnetic flux inside the impenetrable tube 
remains invariable with the variation of the tube radius. 
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Figure 7: The Casimir force as a function of the tube radius. The variable Xq is along 
the abscissa axis, the value of the dimensionless Casimir force ^ } * s along the ordinate 
axis. 

4 Summary 

In the present paper we consider the vacuum polarization effects which are induced in 
scalar matter by imposing a perfectly reflecting (Dirichlet) boundary condition at the 
edge of an impenetrable finite-radius tube which carries magnetic flux lines inside itself. 
Restricting ourselves to a plane which is orthogonal to the tube, we define the induced 
vacuum energy density, see (7) and (8), and analyze numerically its behavior as a function 
of the distance from the tube for the magnetic flux equal to half of the London flux 
quantum ($ = vre -1 ), the tube radius equal to r = 10~ 3 m _1 and different values of the 
coupling to the space-time curvature scalar (£), see Fig. 3. The emergence of the energy 
density, as well as of other components of the energy-momentum tensor, in the vacuum 
can lead to various semiclassical gravitational effects which were estimated under the 
neglect of the tube radius in [15]. 

The present paper summarizes and extends our previous study in [11, 12], and this 
allows us to draw conclusions about the behavior of the total induced vacuum energy, 
i.e. the density integrated over the whole plane, and the Casimir force as functions of the 
tube radius. We find that the total induced vacuum energy is finite and independent of £, 
as long as the tube radius is taken into account. Although the values of the total induced 
vacuum energy are negligible for r ~ m" 1 (see also [12]) being of order 10~ 10 x 27im, 
they are of order 10 _1 x 2nm for r ~ 10~ 3 m _1 , see the Table and Fig.6. The induced 
vacuum energy gives rise to the Casimir force which is directed from the inside to the 
outside of the tube. The force acts at the increase of the tube radius and the decrease of 
the induced vacuum energy under the steady magnetic flux. The force takes considerable 
values at small values of the tube radius and actually disappears otherwise: it is, e.g., 
10 x 2iim 2 at r ~ 10 _3 m _1 and 10 -4 x 27rm 2 at r ~ 10 _1 m _1 . The behavior of the force 
as a function of the tube radius is illustrated by Fig. 7. 

It should be noted that in our case the Casimir force is caused by boundary conditions 
imposed at the boundary enclosing a magnetic flux. The force is periodic in the flux 
value with a period equal to the London flux quantum, attaining its maximal value at 
$ = (2n + l)7re _1 and vanishing at $ = 2mre~ 1 (n G Z). 
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